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Abstract 

Epidemics have so far been mostly studied in undirected networks. However, many real-world networks, 
such as the social network Twitter and the WWW networks, upon which information, emotion or malware 
spreads, are shown to be directed networks, composed of both unidirectional links and bidirectional links. 
We define the directionality £ as the percentage of unidirectional links. The epidemic threshold r c for 
the susceptible-infected-susceptible (SIS) epidemic has been proved to be 1/Ai in directed networks by N- 
intertwined Mean-field Approximation, where Ai, also called as spectral radius, is the largest eigenvalue of 
the adjacency matrix. Here, we propose two algorithms to generate directed networks with a given degree 
distribution, where the directionality £ can be controlled. The effect of £ on the spectral radius Ai, principal 
eigenvector xi, spectral gap (Ai — | A2 1 ) and algebraic connectivity |^ijv-i| is studied. Important findings 
are that the spectral radius Ai decreases with the directionality £, and the spectral gap and the algebraic 
connectivity increase with the directionality £. The extent of the decrease of the spectral radius depends on 
both the degree distribution and the degree-degree correlation po- Hence, the epidemic threshold of directed 
networks is larger than that of undirected networks, and a random walk converges to its steady-state faster 
in directed networks than in undirected networks with degree distribution. 

1 Introduction 

Much effort has been devoted to understand epidemics on networks, mainly because of the increasing threats 
from cybercrime and the expected outbreak of a new fatal virus in populations. Epidemics have been studied 
on undirected networks for a long time and many authors (see [l][2j[3l[4l[5l|6l[7]) addressed the existence of 
an epidemic threshold r c in the susceptible- infected-susceptible (SIS) epidemic process [5]. We consider the 
SIS epidemic process in an undirected network G(iV, L), characterized by a symmetric adjacency matrix A 
consisting of elements that are either one or zero depending on whether node i is connected to j or not. 
Each node i has two possible states at time t: either Xi(t) = meaning that the node is healthy and susceptible 
or Xi (t) = 1 for an infected node. Initially, a certain percentage of nodes is randomly selected as infected. The 
infection from an infected node to each of its healthy neighbors and the curing of an infected node are assumed 
to be independent Poisson processes with rates /? and 5, respectively. Every node i at time t is either infected, 
with probability Vi(t) = Prob[A^(t) = 1] or healthy (but susceptible) with probability 1 — Vi(t). This is the 
general continuous-time description of the simplest type of a SIS epidemic process on a network. The epidemic 
threshold r c separates two different phases of the epidemic process on a network: if the effective infection rate 
t = (3/5 is above the threshold, the infection spreads and becomes persistent; if r < r c , the infection dies 
out exponentially fast. The epidemic threshold ri^ = j- of the SIS model by the N-intertwined Mean-field 
Approximation (NIMFA) [9j [10] lower bounds the real threshold, where Ai is the largest eigenvalue of the 
adjacency matrix A, also called the spectral radius. 
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Table 1: The percentage of unidirectional links in real- world networks 
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Topologies of complex networks have been mostly modeled by Erdos and RenyjJ [XT] as binomial networks 
or by Barabasi and Alber1§] 12: as power-law networks or by Watts and Strogatjf] [13] as small-world networks. 
However, many real- world networks are directed networks as shown in Table [TJ Two kinds of links, namely 
bidirectional links and unidirectional links, exist in directed networks. If node % is connected to node j then 
j is also linked to i, one bidirectional link exists between nodes i and j; and if either i — > j or j — > i exists, 
but not both in between a node pair i and j, a unidirectional link exists. Here, we define the directionality as 
£ = ^unidirectional / f L, where the number of links L = J2i 12 j a ij = u T Au and L = L unidirectiona i+2L bidirec ti ona i, 
because each bidirectional link consists of two unidirectional links with opposite directions. A directed network 
with the directionality £ is denoted by When £ = 0, the network G^ -* is a purely bidirectional network, 

whose adjacency matrix is symmetric, i.e. G^ =0 ) and its corresponding undirected network has the same 
spectral properties. When £ = 1, the network G^ =1 ' is a directed network without any bidirectional link. A 
high directionality is observed in Twitter, as shown in Table [1] A link runs from user A to user B if user 
A follows user B in Twitter, where user A is called the "follower" of user B. Just because user A is user B's 
"follower" does not necessarily mean that the reverse is also true. For example, a famous football player could 
have millions of followers but he may not follow many others. This explains the high directionality £ of Twitter. 
Another type of social networks, the virtual-community social networks, such as Facebook, Chinese Renren and 
LiveJournal, have low directionality. One of the reasons for building this kind of websites is to keep connection 
with real-life friends and classmates. If two people know each other in real-life, they have a high chance to be 
friends with each other, which means that there are more bidirectional links in this kind of networks. In other 
words, the directionality of this type of social networks is low, which also can be observed in Table [TJ The data 
set of the real-world networks is obtained from [14j [XSJ [16] and the description of these networks is attached 
in the Appendix A. Garlaschelli and Loffredo [17] researched the reciprocity in directed networks, where the 
reciprocity is equal to 1 — £. Here, we focus on how the directionality £ influences the epidemic threshold and 
other spectral properties. 

Stimulated by the directed social networks, we study the SIS epidemic threshold in directed networks. 
Percolation theory for directed networks with £ = 1 was firstly developed by Newman and co-workers [18] [19] . 
Boguha and Serrano |20] pointed out that even a small fraction of bidirectional links suffices to percolate the 
network. Moreover, the bidirectional links will speed up some propagation process |171 121] . The epidemic 

1 An Erdos- Renyi random graph can be generated from a set of N nodes by randomly assigning a link with probability p to each 
pair of nodes. 

2 A Barabasi-Albert graph starts with m nodes. At every time step, we add a new node with m links that connect the new 
node to m different nodes already present in the graph. The probability that a new node will be connected to node i in step t is 
proportional to the degree di(t) of that node. This is referred to as preferential attachment. 

3 A Watts-Strogatz small-world graph can be generated from a ring lattice with N nodes and k edges per node, by rewiring each 
link at random with probability p. 
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threshold Tc — XT 01 the SIS epidemic process by NIMFA also holds for directed networks in (52]. This 
motivates us to investigate the influence of the directionality £ on the epidemic threshold = in directed 
networks. Recently, Chen and Olvera-Cravioto [23] proposed an algorithm to generate a directed network with 
given a degree distribution, which is similar to the configuration model. However, Chen and Olvera-Cravioto's 
algorithm in [23. cannot generate a directed network with a precise directionality £. Here we propose two 
algorithms to generate directed networks with a given degree distribution and a given directionality £, where £ 
is a control parameter. 

Topological properties of directed networks, such as the short loops, closure connectivity, degree distributions, 
domination and communities on realistic directed networks have been studied in 24, 25, 26, 27 . Dynamic 
processes on real directed networks have also been researched in [28j [29] . Here, we investigate the effect of 
the directionality £ on the spectral radius Ai, the principal eigenvector x\, the spectral gap (Ai — A2) and 
the algebraic connectivity /Ujv-i in both binomial and power-law directed networks. We call a network a 
binomial @ (or power-law) directed network, if its in-degree and out-degree follow the same binomial (or power- 
law) distribution. Furthermore, we explore the the influence of the linear degree correlation po (also called 
the assortativity) on the epidemic threshold Tc besides the directionality £ in both binomial and power-law 
directed networks. 

Interestingly, the spectral radius Ai of networks G^ =0 ^ is always larger than that of directed networks G^ =1 ) 
when the degree distribution and the assortativity of these networks are the same. It means that the epidemic 
threshold r c in undirected networks is smaller than that in directed networks with the same degree distributed 
and assortativity. Moreover, the decrease of the spectral radius Ai with £ is large when the assortativity is large 
in binomial directed networks, whereas the opposite is observed in power-law directed networks. We conclude 
that a raising of the percentage of unidirectional links increases the epidemic threshold and unidirectional links 
impose an addition constraint on the spread of information/viruses. 

2 Algorithm Description 

The networks mentioned in this paper are simple, without self-loops and multiple links in the same direction. 
One of the main goals in this work is to propose several algorithms for generating a directed network with a given 
directionality £. The in-degree and out-degree preserving rewiring algorithm proposed in this section, is inspired 
by the degree preserving rewiring, which is presented in 30 . Here, we firstly introduce the degree-preserving 
rewiring in undirected networks, that monotonously increases or decreases the linear degree correlation po (also 
called the assortativity), while maintaining the node degrees unchanged. Then, we describe the in-degree and 
out-degree preserving rewiring algorithm and the link resetting algorithm in detail. 

2.1 Degree-preserving rewiring 

The degree-preserving rewiring 130] can either increase or decrease the assortativity of a network: (a) the degree- 
preserving assortative random rewiring: randomly select two links associated with four nodes and then rewire 
the two links such that the two nodes with the highest degree and the two lowest-degree nodes are connected. 
If any of the new links exists before rewiring, discard this step and a new pair of links is randomly selected; (b) 
the degree-preserving disassortative random rewiring: randomly select two links associated with four nodes and 
then rewire the two links such that the highest-degree node and the lowest-degree node are connected, while 
also the remaining two nodes are connected provided the new links do not exist before. Either rewiring step (a) 
or (b) can be repeated to monotonically increase or decrease the assortativity. 

4 For example, an Erdos-Renyi random network is a binomial network with the linear degree correlation (also called the assor- 
tativity) po = 0. A general binomial network could possibly have an assortativity within a large range. 
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2.2 In-degree and out-degree preserving rewiring algorithm (IOPRA) 

We define our in-degree and out-degree preserving rewiring algorithm (termed IOPRA) as follows: randomly 
choose two unidirectional links with four nodes, and rewire the two unidirectional links. In IOPRA, the head of 
one unidirectional link only can rewire with the head of another unidirectional link, in order to maintain both 
the in-degree and out-degree of the four nodes unchanged (see Figure [T|) . If any link exists between the two new 
node pairs before rewiring, or if the directionality of the new network does not increase (or decrease) towards 
the given directionality, discard this step and randomly select another unidirectional link pairs. We illustrate 
the process of IOPRA in Algorithm [T] (see Appendix C). IOPRA actually changes directionality £ of a given 
network G without changing the in and out-degree of each node. It means that the degree distribution of the 
directed network G^' comes from the original network G. If the original network is an undirected network, the 
in-degree sequence is exactly the same as the out-degree sequence in the directed network. 




(a) (b) 



Figure 1: (Color online) The in-degree and out-degree preserving rewiring 

Applying the algorithms IOPRA to an initial network actually randomizes the connections of the original 
network to change the directionality without changing the degree of each node. Hence, if the initial network is 
e.g. an ER random network or a BA power-law network, where the connections are originally random, IOPRA 
changes only the directionality £. However, if we apply IOPRA to e.g. a lattice, the resulting network has not 
only a different directionality but a more randomized structure. 

2.3 Link resetting algorithm (LRA) 

LRA can only generate random directed networks with a binomial degree distribution. We start with an Erdos 
and Renyi network G P (N) with N and p, and use LRA to change the directionality. We randomly choose 
fraction £ of bidirectional link pairs from G P (N). Then, we randomly choose only one unidirectional link from 
each bidirectional link, and relocate the unidirectional links to the places without any link. A directed network 
G^' with directionality £ is obtained (see Algorithm [2] in Appendix C). LRA is used to generate binomial 
directed networks with controllable directionality £. As the links are randomly reallocated, the in-degree and 
out-degree of every node in the new network is not exactly the same as that of the original bidirectional network 
G P (N). However, the directed network obtained by LRA still has the same binomial degree distribution for both 
in and out-degree. In this work, two types of binomial directed networks can be generated: one is generated 
by IOPRA (called the IOPRA binomial directed networks), whose nodes have the same in-degree and out- 
degree; and the other created by LRA (called the LRA binomial directed networks), has different in-degree and 
out-degree sequences, which is more general in real-world networks. 
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( c ) (d) 

Figure 2: Plot of the spectral radius versus the directionality in (a) and (b), as well as the assortativity versus 
the directionality in (c) and (d), in both binomial and power-law directed networks generated by IOPRA. 

3 Spectral properties in directed networks 
3.1 Spectral radius of the directed networks 

The adjacency matrix of the directed network is an asymmetric matrix, whose spectral radius Ai is still real 
by the Perron-Frobenius Theorem (see [31]). We generate the directed networks from the ER (N = 1000, 
p = 2lnN/N) and BA (N — 1000, m = 4) networks gradually using IOPRA, with the directionality £ ranging 
from to 1. The influence of the directionality £ on the spectral radius Ai and the assortativity pu is studied 
in both power-law directed networks and binomial directed networks (see Figure [2J . 

Apart from some wobbles, the spectral radius Ai almost linearly decreases with the directionality £. More- 
over, the assortativity pd of the network fluctuates slightly around 0. We observe that the tiny leaps of spectral 
radius Ai happen when the assortativity po has a rise, which is understandable, because it has been presented 
in [30] that the spectral radius Ai increases with the assortativity pp. Figure [3] examples that the spectral 
radius Ai may increase instead of decrease when the directionality increases due to the assortativity pjj. We 
will study the effect of the assortativity pu on the decrease of the spectral radius Ai with the directionality £ 
in Section @] later. 

With LRA, we generate binomial directed networks with directionality £ from to 1 with step 0.1. The as- 
sortativity pd of all the binomial directed networks generated by LRA is around 0. The effect of the assortativity 
Pd can be ignored here. The spectral radius Ai is calculated in directed networks with different directionality £. 
We performed all the simulations for 10 3 network realizations. The spectral radius Ai is plotted as a function 
of the directionality £ in binomial directed networks with p — 2lnN/N and p = 0.05 in Figure SJ From the 
observation, the spectral radius Ai decreases linearly with the directionality £ with the slop around —1, which 
is independent from the link density p of the networks. We make a proposition which is used to calculate the 
spectral radius Ai of the binomial directed networks with directionality £. 
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4 = 2.7422 4 = 2.7567 

Figure 3: (Color online) Example: the spectral radius increases with the directionality £, because of the increase 
of the assorativity (where pD(Guft) = -0.6190, £(G; e /t) = 0.8333, and p D {G righ t) = -0.5714, £,(G ng ht) = 
0.9167). 
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Figure 4: (Color online) The average of the spectral radius as a function of the directionality for binomial directed 
networks generated by LRA with size N = 1000. Two values for the link density p are shown: p = 2lnN/N 
(circles) and p = 0.05 (diamonds). 

Proposition 1 Let G P (N) be an Erdos-Renyi (ER) random network with N nodes and link density p, and 

be a binomial directed network generated by LRA whose in and out- degree follow the same binomial distribution 

as G P (N). The average of the spectral radius of these two types of networks satisfies 



E[Xi(G^)] ~ E[Xi(G p (N))] - £ 



(1) 



Arguments: 

To analyse the spectral radius Ai of a binomial directed network generated by LRA with link density p, we 
consider that the network is obtained by adding 2p£(^) unidirectional links to a bidirectional ER network with 
size N and link density p(l — £). The average spectral radius of ER networks with N and p is E[Xi(G p (N))] 
: Y — 2)p + 1 + 0{A=). Thus, the average spectral radius of ER network with size N and link density p(l — £) 



is 

EMG^i-qW)] = (N- 2)p(l - + 1 + 0(-!=). 

The principal eigenvector of an adjacency matrix A is denoted by x\ obeying the normalization x\x\ = 1. Let 
C denote the adjacency matrix of the resulting network after adding one unidirectional link to network G. If 
one element in an non-negative matrix A is increased, then the largest eigenvalue also increased [311 pp. 236, 
Lemma 7] as 

Ai(C)> X 1 (A) + (x 1 ) i (x 1 ) j . 

Hence, the average increase of the spectral radius by adding m unidirectional links in random networks is 
obtained as 

£[Ai(C) - Xx(A)} > mElixiUxx^], 
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The sum of the product of components in the principal eigenvector of Erdos and Renyi networks is approximated 
by a function of link density p (see Figure[5]) . When N — > oo, and p ^ 1 or 0, the fitting function can be expressed 
as, 

N N 

(Xi)i{x\)j 

j=l i=l 

Since J2i=ii x i)i = 1 an d using the linearity of the Expectation Operator E[X], we obtain 



= N- -±0(1). 
P 



E [(xx)i(xx)j 



N - i-2 
p 



(2) 



' Ji N(N-l)' 

Hence, the average of spectral radius of the directed network by adding m — 2p^(^) unidirectional links to 
network Gwi_£)(iV) is, 

E[Xi(G {C> )] ~ ^[A 1 (G p(1 _ S) (iV))] + 2 (N(N - l)/2) p£E [(xM^j) • 

Using ©, 



E[Xi(G^)] ~ (N - 2)p(l -0 + 1 + 2 1} ) 
= (7V-2>+l-£, 

which leads to (1). □ 

Juhasz [32] also pointed out that the largest eigenvalue Ai(G^ =1 )) of a random directed network with 
link density p and size N is almost surely Np, when N — > oo. In ER random networks, the spectral radius 
£'[Ai(G^ =0 ')] = (N — 2)p + 1 + 0(^7=) — > iVp + 1, when — > oo. Both earlier results are consistent with our 
proposition and support that the slope of the decrease of the spectral radius Ai with directionality £ is around 
-1. 
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Figure 5: (Color online) Sum of the product of components in the principal eigenvector as a function of the link 
density p in ER networks (TV = 1000). 



With Proposition [T] we analyse the effect of the size iV on the largest decrease of the spectral radius 



A = 



_ A 1 (G^=°))-A 1 (G^ = 1 )) 



Ai(G«=°>) 



We make the prediction that A — ¥ if N — ¥ oo for binomial directed networks, since 



the decrease of the spectral radius is almost a constant value, meanwhile, the spectral radius Ai(G^ =0 - ) ) of dense 
binomial directed network is large and increases with the size of the networks. This implies that the effect of 
directionality £ on the spectral radius is small in large binomial networks. 
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3.2 Principal eigenvector in the directed networks 

In view of the importance of the principal eigenvector x\ in characterizing the influence of nodes on the spectral 
radius, in this section, we explore the principal eigenvector in directed networks. The principal eigenvector of 
the directed networks with the directionality £ from to 1 with step 0.1, are calculated. Then, the components 
in the principal eigenvector are sorted in an ascending order. The simulations for 10 3 network realizations are 
performed. Figure |6] illustrates that the components in the principal eigenvector are more uniform in binomial 
directed networks. The principal eigenvector x\ — » where u is the all-one vector, when the network has 
a higher directionality £. Moreover, the variance of components in the principal eigenvector linearly decreases 
with the directionality £ in both binomial directed networks and the power-law directed networks (see Figur^Z]) . 
The decrease of the variance l^ar[xi] in binomial directed networks by LRA is larger than that in binomial 
directed networks by IOPRA. Thus, nodes in LRA binomial directed networks have more equal contributions 
to the spectral radius than nodes in IOPRA binomial directed networks, when the directionality is the same. 
The difference between LRA and IOPRA binomial directed networks lies in the fact that LRA allows each 
node to have a different in- and out- degree. The in- and out- degree distribution are the same in both LRA 
and IOPRA binomial networks. Hence, the connections in LRA binomial directed networks are more random 
contributing to a smaller yetr[xi]. Li et al. [53] has shown that both a large variance of nodes' degree and a 
large assortativity pu contribute to a large variance Var[xi] of the components in the principal eigenvector x\. 
Here, we point out that the directionality £ is also related to the variance Var[xi] of the components in X\. 
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Figure 6: (Color online) The change of the components of principal eigenvector from bidirectional networks to 
directed networks. 
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Figure 7: (Color online) Plot of the variance of the principal eigenvector versus the directionality (a) in binomial 
directed networks (generated by LRA in square and by IOPRA in triangle) and (b) in power-law directed 
networks by IOPRA (10 3 network realizations). 
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Figure 8: (Color online) Plot of the spectral gap as a function of the directionality (A) in binomial directed 
networks (generated by LRA in square and by IOPRA in triangle) and (B) in power-law directed networks (10 3 
network realizations). 

3.3 Spectral gap of the directed networks 

The difference between the largest eigenvalue Ai and the second largest A2 is called the spectral gap. All 
eigenvalues of the symmetric adjacency matrix of an undirected network are real. Here we focus on the directed 
networks, whose adjacency matrix is asymmetric. The eigenvalues of the directed networks can be complex 
numbers. Normally, the modulus is used to measure the magnitude of an imaginary number. The spectral gap 
Ai — IA2I increases with the directionality £ in both the directed binomial networks and the power- law directed 
networks (see Figured]). The larger the spectral gap is, the faster the random walk converges to its steady-state. 
Thus, the dynamic process in the directed network reaches the steady-state faster than that in the undirected 
network with the same degree distribution. Figured (a) implies that a dynamic process is slightly faster to reach 
the steady-state in IOPRA binomial directed networks than in LRA binomial directed networks. The existence 
of large spectral gap together with a uniform degree distribution results in higher structural sturdiness and 
robustness against node and link failures. Hence, directed networks with high directionality £ and a uniform 
degree distribution are more robust than undirected networks with large variance of degree. 



3.4 Algebraic connectivity of the directed networks 

The Laplacian Q = \BB T = A — A, where the incidence matrix B is an N x L matrix with elements [3T] 

!1 if link ei = i — > j 
— 1 if link e; = j — > i 
otherwise, 

A = |(Aj„ + A out ), A in and A out are the in-degree and out-degree matrices, and A = ^(A + A T ). If the 
network is an undirected network, A is the adjacency matrix A and A = diag(c?i, d,2, ■ ■ ■ >djv) is the degree 
matrix. The second smallest eigenvalue /Ltjv-i of the Laplacian Q was named by Fiedler [34]. Here we study 
the algebraic connectivity in directed networks, The algebraic connectivity is always mentioned together with 
the spectral gap, since they may both be used to quantify the robustness and the network's well-connectedness. 
The larger the algebraic connectivity is, the more difficult it is to cut the network into disconnected parts. Here, 
we study the influence of the directionality £ on the algebraic connectivity //jv-i in directed networks. The 
Laplacian Q of a directed network is asymmetric, thus, the algebraic connectivity hn-i might be a complex 
number and the modulus |//jv-i| is used to measure the magnitude of the algebraic connectivity. We point out 
that, compared to the real part, the imaginary part of the algebraic connectivity is rather small, which can 
even be ignored. The algebraic connectivity increases with the directionality £ in both the binomial directed 
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Figure 9: (Color online) Plot of the algebraic connectivity as a function of the directionality (A) in binomial 
directed networks (generated by LRA in square and by IOPRA in triangle) and (B) in power-law directed 
networks (10 3 network realizations). 



networks and the power-law directed networks (see Figure [9]). It illustrates that the directed networks are more 
difficult to break into parts. Moreover, the LRA binomial directed networks is more strongly connected than 
the IOPRA binomial directed networks with the same in- and out- degree distribution. 



4 Effects of the assortativity on the spectral radius of the directed 
networks 

In Section [3J we have discussed how the spectral properties change with the directionality in directed networks, 
where the assortativity is always around zero. Here, we study further how the spectral radius Ai changes with 
the directionality £ when the assortativity pd is the same, and how the change alters with the assortativity. 
Two approaches are applied to investigate this problem. 

Approach 1: We firstly do degree-preserving rewiring on ER networks (or BA networks) to obtain a set of 
bidirectional networks with assortativity pr> from —0.8 to 0.8 (or —0.3 to 0.3) with step 0.1. Secondly, alter 
the directionality (, of all bidirectional networks with each assortativity using IOPRA. The directionality £ is 
changed from to 1 with step 0.1. IOPRA randomizes network connections, and thus push the assortativity 
of the resulting directed network toward zero if the original network has a non-zero assortativity. Then, the 
spectral radius Ai of all directed networks with £ from to 1 with step 0.1, is scatter plotted as a function of the 
assortativity po in different colorful lines, from the top gray line to the bottom pink line. Figure [TOl plots the 




simulation result of one binomial network realization and 10 2 binomial network realizations. The simulation of 
one realization is almost same as the result of a large number of network realizations, which points to almost 
sure behavior [35] . The results in the power-law directed networks are shown in Figure 1111 
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Figure 11: (Color online) The spectral radius as a function of the assortativity in power-law directed networks. 

Approach 2: Firstly, generate ER networks (or BA networks) G^ - 1 whose directionality £ = 0, then, use 
IOPRA to generate binomial directed networks (or power-law directed networks) G^ =1 ) with directionality £ = 
1. Secondly, do degree-preserving rewiring on ER networks (or BA networks) G^ =0 ) to change the assortativity 
of the undirected networks, and the spectral radius of networks are calculated and plotted as a function of the 
assortativity (squares in Figure [T2|). Then, apply IOPRA to change the assortativity of the binomial directed 
networks (or power-law directed networks) G^ =1 ' . IOPRA used here, takes the assortativity as the judgement 
for whether to do rewiring at each step as long as the directionality remains at £ = 1. The spectral radius Ai of 
the directed networks G^ =1 - ) with different assortativities is also shown as a function of the assortativity (circles 
in Figure [T2|) . 

Observations from Figures [TUl [H] and Q21 all show that the spectral radius Ai always decreases with the 
directionality £ when the networks have the same degree distribution and the same assortativity pr>. Moreover, 
the degree distribution of the network also influences the change range of the spectral radius Ai with £. The 
decrement of the spectral radius Ai with £ increases with the assortativity in binomial directed networks (see 
Figures [TUl and fl2l (a)). On the contrary, the decrement of the spectral radius Ai with £ goes down with the 
assortativity in power-law directed networks (see Figures [Til and [T2l (b)). Furthermore, the decrease of the 
spectral radius in power-law directed networks is larger than that in binomial directed networks, when the 
assortativity is 0. 

Summarizing, the spectral radius Ai decreases with the directionality £ when the assortativity is the same. 
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Figure 12: (Color online) The spectral radius as a function of the assoratativity (a) in binomial directed 
networks (N — 1000, p — 2lnN/N) and (b) in power-law directed networks (JV = 1000, m = 4) for 10 3 network 
realizations. 
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In order to protect the network from the virus spreading via increasing the epidemic threshold, while guar- 
anteeing the degree distribution and the assortativity maintained, increasing the directionality of networks is 
recommended, which, in return enhances the topological robustness. 

5 Conclusions 

In this paper, we illustrate that many real-world networks are directed. Several algorithms used to generate 
the directed networks with given directionality £ are proposed. This allows us to study the influence of the 
directionality £ on the spectral properties of networks. The spectral radius Ai, which is the inverse of the SIS 
NIMFA epidemic threshold rl X \ is studied in directed networks. A universal observation is that, the spectral 
radius decreases with the directionality when the degree distribution and the assortativity of the network is 
preserved. In order to suppress virus spreading throughout a network, increasing the directionality of the network 
is advocated to raise the epidemic threshold. The raising range of the epidemic threshold is also influenced by 
both the degree distribution and the assortativity of the networks. Moreover, the variance of the components of 
the principal eigenvector decreases with the directionality, and the spectral gap and the algebraic connectivity 
increase with the directionality, which illustrates that a rise of the directionality enhances the robustness of 
the networks. We have compared spectral properties in LRA binomial directed networks and IOPRA directed 
binomial networks. The influence of the difference between in- and out- degree of nodes on spectral properties 
for power-law directed networks is an open question. The information dissemination is similar to the virus 
spreading on networks. Our results suggest that encouraging directed online social relations could possibly 
prevent violent emotions (such as in 2011 London riots) from propagation. 
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Appendix 

A. Introduction of real-world networks 

1. Enron: 

This dataset was made public by the Federal Energy Regulatory Commission during its investigations: it 
is a partially anonymised corpus of e-mail messages exchanged by some Enron employees (mostly part of the 
senior management). We turned this dataset into a directed graph, whose nodes represent people and with an 
arc from x to y whenever y was the recipient of (at least) a message sent by x. 

2. Ljournal-2008: 

LiveJournal is a virtual-community social site started in 1999: nodes are users and there is an arc from 
x to y if x registered y among his friends. It is not necessary to ask y permission, so the graph is directed). 
This graph is the snapshot used by Flavio Chierichetti, Ravi Kumar, Silvio Lattanzi, Michael Mitzcnmacher, 
Alessandro Panconesi, and Prabhakar Raghavan in "On compressing social networks", KDD '09: Proceedings 
of the 15th ACM SIGKDD international conference on Knowledge discovery and data mining, pages 219-228, 
2009, ACM press, and was kindly provided by the authors. 

3. Twitter-2010 

Twitter is a website, owned and operated by Twitter Inc., which offers a social networking and microblogging 
service, enabling its users to send and read messages called tweets. Tweets are text-based posts of up to 140 
characters displayed on the user's profile page. 

This is a crawl presented by Haewoon Kwak, Changhyun Lee, Hosung Park, and Sue Moon in "What is 
Twitter, a Social Network or a News Media?" , Proceedings of the 19th International World Wide Web (WWW) 
Conference, pages 591-600, 2010, ACM press. Nodes are users and there is an arc from x to y if y is a follower 
of x. In other words, arcs follow the direction of tweet transmission. 

Note that the distance distribution reported in the paper (Figure 4) is quite different from the one we 
computed (with relative standard error < 0.4% on all points of the cumulative distribution function). Corre- 
spondingly, the average distance reported in the paper (4.12) is quite different from our estimate. We have 
verified that we obtain the distribution reported here even when using breadth-first sampling. 

The node numbering of the original graph was not compact, as the node identifiers were Twitter's actual 
internal identifiers. We thus renumbered nodes in a contiguous way. The original identifiers can be found in the 
file with extension .ids. Thus, you can access the Twitter page associated to a node by getting the corresponding 
identifier and then using the Twitter API. For instance, the node of maximum outdegree is node 2997469, 



corresponding to identifier 19058681, whose owner can be found at http://api.twitter.eom/l/users/show.xml? 
userJd=19058681. 

4. Word Association-2011 

The Free Word Association Norms Network is a directed graph describing the results of an experiment of 
free word association performed by more than 6000 participants in the United States: its nodes correspond to 
words and arcs represent a cue-target pair (the arc x->y means that the word y was output by some of the 
participants based on the stimulus x). 

5. WWW networks 

The networks, "cnr-2000" , "in-2004" , "eu-2005" , "uk-2007-05@100000" and "uk-2007-05@1000000" are small 
WWW networks that were crawled from the Interenet. The "cnr-2000" is crawled from the Italian CNR domain. 
A small crawl of the .in domain performed for the Nagaoka University of Technology is in data "in-2004". The 
"eu-2005" is a small crawl of the .eu domain. This network "uk-2007-05@100000" and "uk-2007-05@1000000" 
have been artificially genereted by combining twelve monthly snapshot of the .uk domain and collected for the 
the DELIS project. 

B. Eigenvalues of the directed networks 
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The spectral radius Ai and the spectral gap (Ai — A2) are considered as important metrics for the percolation 
processes on networks. Here we also give a whole impression of all eigenvalues in directed networks in a Image- 
Real figure. We sort the eigenvalues based on the real part of the eigenvalues on 10 3 simulation realizations. If we 
order them based on the modulus of the eigenvalues, the average of the eigenvalues of 10 3 simulation realizations 
will be balanced. In that case, the results cannot show how the real and image part of the eigenvalues change 
with the directionality £. The changes of the eigenvalues A 2 ; with the directionality £ from to 1 with step 0.1 
in binomial directed networks (N = 10, p = 0.25) are shown on Figure [131 Surprisingly, the real part of all 
eigenvalues tends to 0, when the directionality £ increases. 
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Figure 13: The change of the eigenvalues from bidirected graph to directed graph. 

C. Algorithms 

Algorithm 1 IOPRA(G, £) 

1: Create a bidirectional network G(N, L); 

2: Save network G(N,L) as G s and calculate the directionality 6 of network G s ; 

3: while |6 - f | > 10~ 5 do 

4: Randomly select two unidirectional links i — > j and k — ¥ I associated with the four nodes i, j, k, I; 

5: Rewire the link pair i — > j and k — > I into i — > I and k — > j. The new network G n is obtained; 

6: calculate the directionality of the network G n ; 

7: if |6 - £| > 16, - £| then 

8: G s ^— G n ; 

9 ; 6 ^ £713 

10: else 

11: give up this rewired node pair; 
12: end if 
13: end while 
14: return G s 



Algorithm 2 LRA(G 7 £) 

1: Create a binomial bidirectional network G(N, p); 

2: Randomly choose £ percentage of bidirectional link pairs; 

3: Randomly choose one unidirectional link from each link pair; 

4: Randomly reset the chosen unidirectional links to the locations without any link; 

5: Save the new network as G s ; 

6: return G s 
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